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Time-dependent phase shift of a retrieved pulse in off-resonant EIT-based light storage 
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We report measurements of the time-dependent phases of the leak and retrieved pulses obtained 
in EIT storage experiments with metastable helium vapor at room temperature. In particular, 
we investigate the influence of the optical detuning at two-photon resonance, and provide numer¬ 
ical simulations of the full dynamical Maxwell-Bloch equations, which allow us to account for the 
experimental results. 


I. INTRODUCTION 

Because they do not interact with each other and can 
be guided via optical fibers over long distances with rel¬ 
atively low losses, photons appear as ideal information 
carriers and are therefore put forward as the "flying 
qubits" in most of quantum communication protocols. 
The design of memories able to reliably store and re¬ 
trieve photonic states is, however, still an open problem. 
The most commonly studied protocol, considered to im¬ 
plement such a quantum memory, is electromagnetically 
induced transparency (EIT) [l|. This protocol was imple¬ 
mented in various systems such as cold atoms, gas cells, 
or doped crystals [2|-[4j]. Although the Doppler broaden¬ 
ing might seem to lead to strong limitations, EIT-based 
light storage in warm alkali vapors gives good results and 
is still a subject of active investigation Q. In the last 
years, some experiments were also performed in a Ra¬ 
man configuration, using pulses which are highly detuned 
from the optical resonances in gas cells j6H8|- 

The EIT-based storage protocol in a A atomic sys¬ 
tem relies on the long-lived Raman coherence between 
the two ground states which are optically coupled to the 
excited level. When a strong coupling beam is applied 
along one of the two transitions, a narrow transparency 
window limited by the Raman coherence decay rate is 
opened along the other leg of the system. Because of the 
slow-light effect associated with such a dramatic change 
of the medium absorption properties, a weak probe pulse 
on the second transition is compressed while propagat¬ 
ing through the medium. When this pulse has fully en¬ 
tered the atomic medium, it can be mapped onto the 
Raman coherences which are excited by the two-photon 
process by suddenly switching off the coupling beam. It 
can be safely stored during times smaller than the life¬ 
time of Raman coherence. Finally, the signal pulse can 
be simply retrieved by switching on the coupling beam 
again. In the Raman configuration, the coupling and 
probe pulses are optically far off-resonance but still ful¬ 
fill the two-photon transition condition. The advantage 
is a large bandwidth, that allows to work with data rates 
higher than in the usual EIT regime [6]. 

Atoms at room temperature in a gas cell are partic¬ 
ularly attractive for light storage because of the sim¬ 
plicity of their implementation. The effects of the sig¬ 
nificant Doppler broadening can be minimized using co¬ 
propagating coupling and probe beams, so that the two- 


photon resonance condition can be verified for all velocity 
classes: all the atoms can thus participate to the EIT phe¬ 
nomenon as soon as they are pumped in the probed level. 
As a consequence, handy simple gas cells have turned out 
to be attractive for slow or even stopped light experi¬ 
ments Q. In a previous work Q, we have reported on an 
added phase shift recorded for EIT-based light storage 
experiments carried out in a helium gas at room tem¬ 
perature when the coupling beam is detuned from the 
center of the Doppler line. The simple model that we 
have derived could not satisfactorily account for our ob¬ 
servations that were recorded for intermediate detunings, 
e.g. close to the Doppler broadening of the transition. In 
the present paper, we come back to this problem and pro¬ 
vide new experimental results, i.e. time-dependent mea¬ 
surements of the retrieved signal phase shift, as well as 
numerical results obtained through the simulation of the 
full system of Maxwell-Bloch equations. The behaviour 
of these phase shifts with the coupling detuning seems 
satisfactorily accounted for by our simulations. We also 
perform numerical calculations in the Raman regime. 

The paper is organized as follows. In Section [XT] we 
present the system and setup and describe how to mea¬ 
sure the time-dependent phase shift of the retrieved pulse 
with respect to the coupling beam. We also briefly re¬ 
call the system of Maxwell-Bloch equations which gov¬ 
erns our system and describe their numerical integration. 
In Section uni we provide our experimental and numer¬ 
ical results and show that they qualitatively agree. We 
also apply our simulations to the far off-resonant Raman 
case. Finally, we conclude in Section HVl and give possible 
perspectives of our work. 


II. EXPERIMENTAL SETUP AND 
NUMERICAL SIMULATIONS 

A. EIT storage experimental setup 

The atoms preferably used for EIT storage experiments 
are alkali atoms, mainly rubidium and sometimes sodium 
or caesium. We choose here to work with metastable 4 He 
atoms, which have the advantage of a very simple struc¬ 
ture without hyperfine levels: transitions are thus far 
enough one from another to investigate the effect of de¬ 
tunings of the coupling and probe beams on light storage 
and retrieval. 
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In our setup represented in Fig. [Q a 6-cm-long cell is 
filled up with 1 Torr of helium atoms which are con¬ 
tinuously excited to their metastable state 2 3 Si by a 
radio-frequency (rf) discharge at 27 MHz. Each of the 
metastable ground states |2 3 S'i, mj = 0, ±l) is hence fed 
with the same rate, denoted by The cell is isolated 
from magnetic field gradients by a three-layer p,-metal 
shield to avoid spurious dephasing effects on the different 
Zeeman components. A strong circularly-polarized field, 
called the control beam, propagates along the quantiza¬ 
tion axis 2 . Its power is set at 18 mW for a beam diameter 
of 3 mm. As shown in Fig. [2J the coupling field drives 
the transitions |2 3 S’i, mj = —l) «-»■ |2 3 Pi , mj = 0) and 
|2 3 5i,toj = 0) |2 3 Pi,mj = l). Owing to the spon¬ 
taneous transitions |2 3 Pi,toj = 0) —> |2 3 5i,toj = ±l) 
and |2 3 Pi,toj = l) —> |2 3 Si,toj = 0, l), the atoms end 
up in the state |1) = |2 3 Si, mj = l) within a few pump¬ 
ing cycles after the coupling beam has been switched 
on. As the atoms are at room temperature, the Doppler 
broadening in the cell is Wd/2tt « 1 GHz. We denote by 
A c = uj c — cuq the detuning of the coupling frequency 
lu c with respect to the natural frequency wo/( 27 r) « 
2.8 x 10 14 Hz of the transition 2 3 S' 1 <->• 2 3 P 1; at the center 
of the Doppler line. 

Once optical pumping is achieved, a weak signal pulse 
is sent through the atomic medium along the z axis. Its 
polarization is circular and orthogonal to that of the cou¬ 
pling beam: the signal therefore couples the state |1) to 
|e) = |2 3 Pi, mj = 0) and we denote by A s = w s —wo the 
detuning of the signal frequency u> s from the center of the 
Doppler profile. Both signal and coupling beams are de¬ 
rived from the same laser diode, and their frequencies and 
amplitudes are controlled by two acousto-optic modula¬ 
tors. Due to the efficiency of optical pumping through the 
coupling beam, we assume that the state |2 3 5 'i, mj = 0) 
remains essentially unpopulated during the whole process 
and we accordingly neglect the driving of the transition 
|2 3 Si, toj = 0) O |2 3 Pi ,mj = —l) by the signal field. 
Submitted to the coupling and signal fields, the atoms 
therefore essentially evolve in the three-level A system 
{|—1} = |2 3 Pi, mj = -l) , |e), |1)} (see Fig. 0 as long 
as the detunings A CjS = w CiS —wo, of the coupling and sig¬ 
nal fields respectively, are small enough to avoid exciting 
neighbouring transitions. Thanks to the absence of hy- 
perfine structure, the range of allowed values for A c s is, 
however, much larger than in alkali vapor experiments: 
indeed, on the positive detuning side the nearest state 
( 3 P 0 ) is 30 GHz away from optical resonance A Ci „ = 0, 
while, on the negative detuning side, the nearest state 
( 3 P 2 ) is 2.29 GHz away. 


Under EIT conditions, the coupling beam opens a 
transparency window for the weak signal beam, which 
can therefore propagate without absorption through the 
medium if its spectrum is not too wide El- In the ex¬ 
perimental results we present hereafter, we used a signal 



Figure 1: Experimental setup for EIT storage in metastable 
helium. The coupling and signal beams are derived from the 
same laser diode. They are initially linearly and orthogo¬ 
nally polarized, of optical frequencies u) c and ui s and Rabi 
frequencies S7 C and fi s , respectively. Acousto-optic modula¬ 
tors are used to control the frequencies and amplitudes of the 
beams. Polarizing Beam Splitters (PBS) allow to separate 
or recombine the beams. Circular orthogonal polarizations 
are obtained by a quarter wave-plate. The cell is contained 
into a /x-metal shielding. After the cell, polarization optics 
select mainly the probe beam with some remaining coupling 
beam.The phase of the coupling beam is scanned thanks to a 
mirror on a piezo-electric tranducer. 


pulse, which consists of a smoothly increasing exponen¬ 
tial followed by an abruptly decaying exponential of re¬ 
spective characteristic times 2 /xs and 150 ns. Its maxi¬ 
mum power is about 170 pW and the beam diameter is 
about 3 mm. Different dissipative mechanisms influence 
the width of the EIT window besides spontaneous emis¬ 
sion, such as collisions and transit of the atoms in and 
out of the beams. These phenomena result in the de¬ 
cay of atomic coherences at the rates 7_n/27r ss 14 kHz 
for the Raman coherence er_n and 7 e i/27r « 22.8 MHz 
for the optical coherence cr e i. We have shown previously 
that velocity changing collisions redistribute the pump¬ 
ing of the atoms over an effective width sligthly smaller 
than the Doppler linewidth (l3 |. In our conditions, with 
a coupling power of 18 mW, this effective width is ex¬ 
perimentally estimated to be Tu/^t: ~ 0.8 GHz. Due to 
power broadening, the width of the transparency window 
is then of the order of 500 kHz. 

The highly dispersive character of the medium under 
EIT conditions can be used furthermore to store and re¬ 
trieve a weak signal pulse: due to EIT dispersion, the 
signal pulse indeed travels with a reduced group velocity 
and is temporally contracted. Once the pulse has entered 
the cell, the coupling beam can be switched off: informa¬ 
tion about the signal pulse is then stored into the Raman 
coherences. After a storage time T, the control beam 
is switched on again, which releases the signal from the 
atomic coherence and ensures EIT absorption-free propa¬ 
gation. Note that in our experimental setup, the optical 
depth is only about 3.5 and the pulse can thus not be 
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Figure 2: Atomic structure scheme for the D1 transition of 
metastable helium. The relevant states which constitute the 
three-level A system are shown in black. A c and A s are the 
optical detunings, and fl c and f i a the Rabi frequencies of the 
coupling and signal beams, respectively. The two photon res¬ 
onance is achieved for A c = A s . 

fully compressed in the cell. Due to the finite width of 
the transparency window and the finite length of the cell, 
a part, typically 10 % of the incoming signal energy, leaks 
out before the coupling beam is turned off and the stor¬ 
age period begins Q. A typical experimental record is 
given in Fig. [3):>. The first part of the detected signal 
is the leak and after a storage time T ss 0.6 /xs, once the 
coupling beam is turned on again, the retrieved signal is 
clearly visible. 

B. Phase measurement setup 

In Q, we investigated the relative phase A <f> of the sig¬ 
nal with respect to the coupling beam, and showed the 
existence of an optical detuning dependent extra phase 
shift ipeit between the incident and retrieved pulse. This 
quantity can be measured through mixing the signal 
emerging from the cell with a small fraction of the con¬ 
trol beam, via polarization optics. The resulting intensity 
thus takes the form 

I = I c + Is + cty/ I C I S (t) cos [A0]. (1) 

In Eq. (fTj). the contrast factor a, which ideally equals 
2, accounts for non-perfect alignment of the beams and 
is measured for each set of data. I c denotes the inten¬ 
sity of the small fraction of the coupling field which is 
mixed and interferes with the signal field. It takes the 
same constant value I c during the writing and retrieval 
periods, while it vanishes during the storage time. The 
value I c is measured in the absence of the signal (one as¬ 
sumes that the introduction of the signal pulse does not 
substantially affect the measurement of J c ). The phase of 
the coupling beam is varied via a piezoelectric actuator 
from one experimental run to another: the scan is slow 
enough so that it is assumed to be constant during both 
the writing or retrieval steps. I s (t) is the time-dependent 
intensity of the signal beam emerging from the cell. We 
denote by /,D (t) and Is ^ (t) the intensities of the leak and 



Time (ps) 


Figure 3: (Color online) Experimental plots showing the data 
processing steps. The coupling beam is switched off between 
times t = 0 and t = 0.6 (is. The optical detunings are set at 
A c ,s = 1 GHz. The same vertical scale is used in graphs a, b 
and c. a) Accumulated plot of the interference signal I (t) be¬ 
tween the probe and coupling beams (blue). Upper and lower 
envelopes I± (t) are shown in orange and red, respectively, b) 
Signal intensity I s (t) at the exit of the cell, deduced from the 
previous accumulated plot a) and from a measurement of the 
coupling intensity I c . c) Interference signal J(t) between the 
probe and coupling beams (blue). It is contained between the 
upper and lower envelopes I± ( t ), shown in orange and red, 
respectively. One can note the presence of spurious oscilla¬ 
tions generated by the acousto-optic modulators, d) Relative 
phases A (black) and A rj>^ r \t) (purple) between the 
signal and coupling beams on the writing and retrieval parts, 
respectively. Notice that A (jr’it) is indeed constant. 


retrieved pulses, respectively. Accordingly, we introduce 
A <jyD (t) and A</>( r ) (f), the relative phases of the leak and 
retrieved pulses, with respect to the coupling beam. 

To obtain the extra phase shift <Peit between the 
incident and retrieved pulses, we measure the relative 
phases AAby homodyne detection. Repeating 
the same writing-storage-retrieval sequence for many dif¬ 
ferent postions of the piezoelectric actuator, we obtain 
an accumulated plot whose upper/lower envelopes cor¬ 
respond respectively to I± (f) = I c + I s (t) ± a^/I c I s (t) 
(see Fig. [3Jr). Given the previously measured value of I c , 
one can infer I s (t) from (/+ + /_) and a from (/+ — J_) 
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(see Fig. GJj)- For a given position of the piezo-actuator, 
one can then obtain Acf)^ ( t ) and A (j)^ (t) through ht- 
ting the experimental record with Eq. 0 at each time 
t (see Figs. d). It was verified both experimentally 
and numerically that the phase of the leak is constant 
(A (f>^ l \t) = Aand can therefore be taken as a ref¬ 
erence for the time-dependent relative phase of the re¬ 
trieved pulse A <j>( r \t). This time independence of the 
leak phase is ensured by the fact that its spectral con¬ 
tent is much narrower than the EIT bandwidth. This 
is obtained thanks to the shape of the signal pulse: a 
slow exponential increase, followed by a sharp decrease. 
The part of the pulse which contains only low frequen¬ 
cies enters first and gives form to a leak, whose phase 
is constant. The “extra phase shift” is then measured 
as [A</>( r ) (f) — A</>(^] = (fiEirit ). Let us stress that in 
0 , we assumed that the relative phases Aof the 
leak and retrieved pulses were time-independent: there¬ 
fore, we directly extracted effective “averaged” values for 
a and 1 -Peit by performing a two-parameter fit of the data 
with Eq. ©■ Here, by contrast, we measure the time- 
dependence of the phases and provide experimental plots 
for ipEirit), without any assumption on its behaviour. 


C. Numerical simulation principles 

For numerical simulations, we described the system in 
the one-dimensional approximation. On the dimensions 
of the atomic sample, the coupling and probe transverse 
profiles are assumed to remain constant. These fields can 

I 


therefore be cast under the form 


E c ,s(z,t) = £ c ,s(z,t)e±e l( “ 


S t fcc 


where £ c s ( z,t ) denote the respective slowly-varying am¬ 
plitudes of the control and signal fields, w CjS and k CtS 
stand for their respective frequencies and wavenumbers, 
while e± = (e^ ± ie y )/-\/2 (e x and e y define an arbitrary 
basis in the plane perpendicular to the propagation di¬ 
rection e^). 

Following e.g. 0 , we model the atomic sample as 
a continuous medium of uniform linear density n a t, and 
define the average density matrix of the slice [z, z + A z] 
by 


a{z,t) 


1 

riatAz 


Z w ■ 

z<Zi<z-j-Az 


We moreover define the density matrix elements (Ty = 
(i \d\ j) where i, j refer to the atomic levels and can take 
the values —1,1 or e (see Fig. (2). We introduce the 
slowly-varying coherences <r e i, tr e _i and cr_n defined by: 


~ iu; s fi—— ) 

cr e i — e V c 'cr e i, 

_ Aujcft— — 1 ^ 

<J e —1 = e ' 

_ iA R(t-z) 
cr-ii = e t c 'er_n, 

and write the Bloch equations in the rotating wave ap¬ 
proximation, for the class of velocity which is at the cen¬ 
ter of the Doppler profile: 


dtcr- 1-1 
(); (lee 
<9 t CTii 

<9*d_n 
<9*0-_le 
<9* die 


dr 1 bo . //-j* ~ 1 1 - \ 

7t^—1 — 1 1 ~^~^ee 1 1 \\l c G e — i \L C (J— le) : 


(r 0 + '"/t) &ee H - i {il c O-le 1 “1“ 

2 7t , r 0 . /y~»* ~ a ~ \ 

.3 7* (7 11 ' ^ le) 5 




cr e 1) 


- ( 7-11 - iA fl ) cr_n + i (D*d e i - D s d_ le ), 

- (7_ie + iA c ) d_ie + i [D* (cr ee - ct_i_i) - n*d_ii], 

- (7ie +iA s ) die + i {£9* (o-ee “ <7ll) ~ . 


( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 
(7) 


Here, To is the population decay rate of the state |e), 
and the Rabi frequencies H C)6 are defined by 

hfl c ,s = 7^dc,s^*c,S ('S'; t) 5 


where d C:S 


d.e± 


T 1 


are the relevant matrix ele¬ 


Doppler width T d • This gives satisfactory results thanks 
to the redistribution of the pumping by velocity chang¬ 
ing collisions fl5i. flQ| . All our simulations were performed 
using this purely homogeneous broadening model. Con¬ 
sequently, we call A C)S the optical detuning, implicitly 
defined with respect to the center of the Doppler profile. 


ments of the dipole operator d. 

To take into account all the atoms that are distributed 
in different velocity classes over the Doppler linewidth, 
we developed a simple model, in which the optical coher¬ 
ence decay rates 71 e = 7 _i e are replaced by the effective 


To ensure that the full population remains constant, 
the discharge-assisted ground-state feeding rate A has 
been set to 7 *, the transit rate of the atoms through the 
laser beam. Moreover, while the state | —1) is fed with the 
rate j = the state 11 } is effectively fed with the rate 
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~ = r] 1 - As can be checked by considering the full 

6-level system including not only the A system of interest 
but also the states |2 3 Si,mj = 0) , |2 3 Pi, mj = ±l), the 
state 11) is indeed directly fed by the rf discharge with the 
rate but also indirectly via the state \2 3 Si,mj = 0) 
whose population is (almost) immediately transferred to 
|1) through optical pumping. 

Finally, in the medium, the fields propagate according 
to the Helmholtz equation, written in the slowly-varying 
envelope approximation 

( d 1 d \ 

V ~dz + cdt ) c,s ( Z,t ' = lrtc ’ s(Te * 1 > ( 8 ) 

where r] C:S = (n at k C:S \d CtS \ 2 ) / (2he 0 ) ■ 

The set of Maxwell-Bloch equations Eqs J2][8]) was nu¬ 
merically solved in Matlab using the Lax discretization 
method 0- The medium was split into 100 spatial steps 
of 0.6 mm while the whole storage/retrieval sequence was 
split into 6 x 10 6 timesteps of 2 ps. We present and dis¬ 
cuss our numerical results in the following section. 

III. RESULTS AND DISCUSSION 

All the experimental results and simulations are per¬ 
formed at two-photon resonance, which means that the 
coupling and signal optical detunings are equal. Fig. 3] 
shows experimental records for the time-dependent extra 
phase shift <pEiT(t), achieved with different values of the 
optical detunings A = A c = A s , between 0 and 2 GHz. 
The detunings are set here on the positive side, where the 
nearest state ( 3 Fo) is nearly 30 GHz away. Each curve is 
obtained after averaging over 15 sets of data, recorded at 
different times, for different positions of the homodyne 
detection piezo-actuator. The traces recorded on the os¬ 
cilloscope present some spurious oscillations at a period 
of about 90 ns. This noise is generated by the acousto¬ 
optic modulators and could be removed by numerically 
filtering the spurious frequencies during the data process¬ 
ing. In Fig. 31 the time origin corresponds to the begin- 
ing of the retrieval, when the coupling beam is turned on 
again. At that time, the probe intensity starts increasing 
to form the retrieved pulse. We only plot the evolution of 
the extra phase shift < Peit{1 ) when the signal intensity is 
high enough, typically from roughly 100 ns to 1 fis after 
the start of the retrieval. One can see that fEirit) is not 
constant over the retrieval, and its magnitude increases 
with the optical detuning A. 

The experimental plots are compared with numerical 
simulations of the full Maxwell-Bloch equations derived 
as explained in the previous section. These simulations 
are in good agreement with experimental results: both 
present the same general shape for tpEiT (i), and the same 
qualitative behaviour with the optical detuning A. One 
possible source for the observed discrepancies is our over¬ 
simplified treatment of the velocity distribution in the 
Doppler profile. Here, we indeed assumed that velocity 
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Figure 4: Temporal evolution of the extra phase shift tfiEiT{t). 
Experimental data, recorded from bottom to top: for A = 0.2 
GHz (blue, bottom), 0.6 GHz (red), 1 GHz (green), 1.3 GHz 
(black), 1.7 GHz (cyan, top). Corresponding simulations are 
shown as continuous lines. The time origin corresponds to 
the begining of the retrieval, when the coupling is switched 
on again. The storage time is set at T = 0.6 ps. 


changing collisions are efficient enough to instantaneously 
and perfectly redistribute atoms pumped in the probed 
level |1) over the effective Doppler profile so that all these 
atoms contribute coherently to the storage process as if 
the broadening were homogeneous. Although this ap¬ 
proximation is commonly used (see 00 ), it might be 
severely questioned here, especially in optically detuned 
conditions that change the thermal equilibrium. In par¬ 
ticular, the absorption of the coupling beam measured 
experimentally could not be well reproduced by the sim¬ 
ulations. This should also have an effect on the storage 
efficiency and on the temporal shape of the phase. Note 
that in the simulation program, in order to “minimize” 
this problem, we use an averaged coupling intensity over 
the length of the cell as the input parameter, instead of 
the real coupling intensity measured at the entrance of 
the cell. Note also that we have checked our numeri¬ 
cal results agree with the analytic approximate solutions 
presented in 0. 

To understand better the physical origin of the extra 
phase-shift ifiEirit), we compared the time-dependent 
relative phase A(f> (t) between the signal and coupling 
beams at the exit of the cell, obtained: i) during the stor¬ 
age and retrieval of a weak signal pulse (A = A <f>^ (t) 
before t = 0 and A <j){t) = A 4^ r \t) after the storage time 
T = 0.6 /is), and ii) during the direct EIT-propagation 
of the same weak pulse in the medium, while the cou¬ 
pling amplitude remains constant. To compute A <p(t) in 
case i), we used the full simulation of Maxwell-Bloch set 
of equations, whereas in case ii) we simply propagated 
each spectral component ui of the incoming pulse with the 
corresponding susceptibility \ (w). Fig. [^simultaneously 
displays the results we obtained in both cases, for two dif- 
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Figure 5: Temporal evolution of the relative phase A cf)(t) dur¬ 
ing the direct propagation of a weak signal pulse through 
the medium under EIT conditions (“Simple propagation”) and 
during the storage and retrieval of the same pulse (“Storage 
sequence”), for two different optical detunings: a) 1.7 GHz 
and b) 0.2 GHz. The time origin is arbitrarily chosen as the 
starting of the storage period. 

ferent values of the optical detuning A = 0.2, 1.7 GHz. 
The shape and order of magnitude for A 4> ( t ) are clearly 
the same in cases i, ii): the main effect of the storage is 
to introduce a delay corresponding to the storage time 
T . Here we chose T = 0.6 [is, but we checked both ex¬ 
perimentally and theoretically that this phase shift does 
not depend on this storage time. This suggests that the 
observed extra phase shift <peit is essentially due to the 
propagation under EIT conditions. Indeed, the stored 
part presents a sharp decrease associated with many fre¬ 
quency components, and is thus very sensitive to disper¬ 
sive effects. This problem should be taken into account 
for high speed information applications, like experiments 
performed in the Raman configuration [ti|-@|. We have 
thus performed simulations in the far detuned regime as 
presented on Fig. [6] The simulation results shown here 
are plotted for three different optical detunings A = 10 
GHz, A = 15 GHz, A = 20 GHz, much higher than the 
T£> « 0.8 GHz Doppler broadening. They were obtained 
for a number of atoms n a t which is 10 times higher than 
in our experimental case, and for a coupling power of 
200 rnW. These simulation results demonstrate a similar 
effect on the retrieved signal pulse phase, even slightly 
stronger than in our experimental conditions. We have 
also checked that our numerical results in the Raman con¬ 
figuration agree with the analytic approximate solutions 
presented in [10]. 

IV. CONCLUSION 

In this paper, we have experimentally investigated 
a time-dependent extra phase shift that appears in a 


storage-retrieval experiment, performed in a room tem¬ 
perature atomic cell, in optically detuned conditions. 



Figure 6: Temporal evolution of the phase shift ifiEirit), for 
different optical detunings in the Raman regime: A = 10 GHz 
(blue, full line), A = 15 GHz (red, dashed line), A = 20 GHz 
(green, dotted line). These curves are obtained with the full 
simulation of Maxwell-Bloch set of equations, for ten times 
more atoms than in our experimental case, and for a coupling 
power of 200 rnW. The time origin is chosen as the starting 
of the retrieval period. 

This phase shift varies with time and does not depend 
on the storage time. We have provided numerical simu¬ 
lations which qualitatively agree with the experimental 
results: in particular, it appears that the magnitude of 
the relative phase depends on the optical detuning, while 
its temporal shape is mainly given by the spectrum of the 
incoming pulse. We explain the existence of this extra¬ 
phase by propagation effects that can be understood by 
a simple propagation model under EIT conditions with 
optically detuned beams. Discrepancies may be due to 
an approximate treatment of Doppler broadening in the 
cell. 

The results presented here might be of importance, 
particularly for light storage experiments performed in 
the far-detuned Raman regime, as reported in [7]. How 
these results translate into the regime of quantum light 
is an intriguing feature that we intend to address in a 
future work. 
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